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' Abstract. We show that the operatorial framework developed by Voiculescu for free 

random variables can be extended to arrays of random variables whose multiplication 
Qj ' imitates matricial multiplication. The associated notion of independence, called ma- 

i tricial freeness, can be viewed as a generalization of both freeness and monotone 

independence. At the same time, the sums of matricially free random variables, called 
random pseudomatrices, are closely related to Gaussian random matrices. The main 
results presented in this paper concern the standard and tracial central limit theorems 
' for random pseudomatrices and the corresponding Hmit distributions which can be 

. viewed as matricial generalizations of semicirle laws. 

> 

c3 . 1. Introduction 

It has been shown by Voiculescu [22] that free random variables arise naturally as 
limits of random matrices. In particular, if we take symmetric matrices whose entries 
> [ form a family of independent Gaussian random variables and we let the size of these 
matrices go to infinity, their moments (with respect to normalized trace composed with 
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OO 



■rj" ■ classical expectation) converge to the moments of freely independent random variables 
, with the semicirle distribution obtained by Wigner [25] as the limit distribution of one 
^ ! Gaussian random matrix. 

00 ! Therefore, we can study free random variables, using at least two different frame- 

O ■ works: operator algebras and random matrices. However, it is to some extent surpris- 
' iiig that a random matrix framework, of quite different nature than that of operator 
^ . algebras, exists for free random variables, and the connection between these two ap- 
proaches does not seem to be very transparent. In this connection, our first motivation 
is to better understand the relation between the operatorial approach to free probability 
and random matrices. 

The second motivation comes from the question whether different types of indepen- 
dence, like freeness and monotone independence of Muraki [16], can be included in 
one natural model. Note in this context that models with more than one state on a 
given algebra, like conditional freeness of Bozejko and Speicher [6] and freeness with 
infinitely many states of Cabanal-Duivillard and lonesco [9,10] extend free probability 
and include, as shown by Franz [11], certain elements of monotone probability. For in- 
stance, this can be done for convolutions, but including monotone independence in the 
framework of conditional freeness can be done only under additional (rather restrictive) 
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assumptions on the considered algebras. We show in this paper that one can remedy 
this situation by introducing a concept of 'independence' which reminds freeness, but at 
the same time has some 'matricial' features which places it somewhere between freeness 
and the model of random matrices. 

A different reason to look for a new concept of independence arises from concrete 
examples of interpolations between free probability and monotone probability [14,15]. 
In particular, the continuous (p, g)-Brownian motions with Kesten distributions and 
related Poisson processes studied in [2,15] lead to the first example of a two-mode 
interacting Fock space introduced directly and not by means of orthogonal polynomi- 
als. This example has certain 'matricial' features which also call for a new model of 
'independence' that would be related to freeness. 

The main result of our paper is the construction of a model, called matricial freeness, 
which is related to the concept of the free product of states introduced and studied by 
Ching [10] in the context of von Neumann algebras and by Avitzour [3] and Voiculescu 
[20] in the context of C*-algebras. The underlying concept is that of the matricially free 
product of an array of Hilbert spaces with distinguished unit vectors {Hi,j,ii,j) which 
reminds the free product of Hilbert spaces. However, it also satisfies the condition 
imitating matrix multiplication and the condition of 'diagonal subordination' which 
says that tensor products must end with 'diagonal' Hilbert spaces. Similarities between 
free probability and 'matricially free probability' hold also on other levels, some of which 
are studied in this paper. Strictly speaking, however, one needs to take a restriction 
of the matricially free product of states, called the strongly matricially free product of 
states, to recover the free product as a special case. 

Roughly speaking, if is an array of self-adjoint matricially free random 

variables in a *-algebra An, equipped with a distinguished state 0„ and a sequence 
(^nj)i^i^n of additional states, called 'conditions', for each natural n, we can take the 
sums 

n 

S{n) = 

called random pseudomatrices, and study their asymptotic distributions (as n ^ co) 
with respect to the states 0„ and with respect to normalized traces 

1 " 

respectively. We assume that the distributions of the Xij(n) in the states 0„ and iljn,k 
are not identical and depend on n in such a way that their ^^-distributions remind 
those of Gaussian random matrices in the approach of Voiculescu. This relation to 
random matrices as well as asymptotic matricial freeness of random pseudomatrices 
will be studied in a forthcoming paper. In progress is the work on the addition of 
matricially free random variables and the resulting matricial R-transform. 

It turns out that the limit theorem for the above random pseudomatrices in the 
states (pn, which we call 'standard', may be viewed as an analog of the central limit 
theorem for free random variables (especially, if we take square arrays). In turn, the 
limit theorem for random pseudomatrices in the states V^n, called 'tracial', is related 
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to the limit theorem for random matrices (especially, if we take square arrays). The 
limit distributions play then the role of multivariate generalizations of the semicircle 
distributions. Let us point out, however, that when we consider triangular arrays, the 
framework of matricial freeness can as well be viewed as a generalization of that of 
monotone independence. 

In Section 2, we introduce the concepts of the 'matricially free product of states' and 
the 'matricially free Fock space'. We obtain from these structures their strong coun- 
terparts in Section 3. In Section 4, we introduce the notions of 'matricial freeness' and 
'strong matricial freeness' and discuss the example of the discrete (strongly) matricially 
free Fock space. In Section 5, of combinatorial nature, we define and study certain 
real-valued functions on the set of non-crossing partitions, defined in terms of traces of 
certain matrices. In Section 6, we study the asymptotic behavior of random pseudoma- 
trices and we prove standard and tracial central limit theorems. The limit distributions, 
which can be interpreted as matricial multivariate generalizations of semicirle laws, are 
studied in Section 7. Their decompositions in terms of s-free additive convolutions in 
the case of two-dimensional arrays are proved in Section 8. Two geometric realizations 
of the limit distributions, in terms of walks on weighted binary trees and in terms of 
weighted Catalan paths, are given in Section 9. 

2. Matricially free products 

In this Section we introduce the notion of the matricially free product of states as 
well as the corresponding notions of the matricially free product of Hilbert spaces and 
the matricially free Fock space. 

When speaking of arrays indexed by two indices, say we shall usually assume that 
i,j e I, where / is an index set. This refers to the situation when we deal with square 
arrays. However, we also want to consider other arrays, like triangular arrays of the 
form T := {{k,l) : k ^ l,:k,l e I}, where / is a linearly ordered index set. Therefore, by 
an array we will understand a subarray of a square array which includes the diagonal. 
Without loss of generality we can use the square array formulation most of the time. 
Of special interest will be the finite-dimensional case when / = [n] := {1,2,..., n}. 

Definition 2.1. Let H := (Tij^) be an array of complex Hilbert spaces. By the 
matricially free Fock space over Ti we understand the Hilbert space direct sum 

m=l (ii,i2)5'---#(jm,im) 
ni,...,nmel^ 

where Q is a unit vector, with the canonical inner product. 

Let us observe that the Hilbert space tensor powers which appear in the above direct 
sum have the following three properties: 

(1) the 'matricial property' - the second index of the preceding power agrees with 
the first index of the following power, 

(2) the 'freeness property' - the consecutive pairs of indices are different, 

(3) the 'diagonal subordination property' - the last pair is 'diagonal'. 
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Of course, if the index set / consists of one element and thus li. is just one Hilbert 
space H., the corresponding matricially free Fock space reduces to the usual free Fock 
space TiTi). In general, however, M-iTi) is a (usually, proper) subspace of the free Fock 

space = '-ij^CHij)- 

Related to the 'matricially free Fock space' is the 'matricially free product of Hilbert 
spaces'. The terminology parallels that introduced in free probability [20,23]. 

Definition 2.2. Let {'Hij,Ci,j) be an array of Hilbert spaces with distinguished unit 
vectors. By the matricially free product of {"Hijj^ij) we understand the pair (Ti.,^), 
where 

00 

with j = Hij QC^ij and ^ being a unit vector. We denote it (H, ^) = 
Proposition 2.1. It holds that 

Proof. We use the definition of the matricially free Fock space, the isomorphism 
A4{7iij) = J-{Hij) for any i,j and regroup terms. ■ 

For any j, introduce diagonal subspaces of H of the form 

CO 

n{j, J ) = ® @ © nl^ (8) (8) . . . (8) 7i^^,^ , 

m=2 (j,i2)i=.-.i=(im,im.) 

and the associated diagonal partial isometrics Vjj : Hjj ®'H{j,j) T~L: 

JiJ ^ J^Jl Jm^Jm^ JiJ J^Jl JmiJm 

where m > 1. 

For any i ^ j we introduce non-diagonal subspaces of Ti. of the form 

CO 

m 3) = ® ® nl, ® (8) ... (8) nl,^ 

and the associated non-diagonal partial isometrics Vij : Hij ®'H{i,j) -^Ti. ioi i j: 
6. (8) (%°,, (8)... (8)^1,- J - %-..®---®^L.v. 

where m ^ 1. 

Each T-t{i,j) is spanned by simple tensors which do not begin with vectors from 7i°j 
and for that reason it is suitable for the left free action of the operators creating such 
vectors. Thus, roughly speaking, both types of partial isometrics jointly replace the 
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unitary maps used in free probability. It is the diagonal subordination property which 
is responsible for distinguishing two types of isometries. 

Consider an array of C*-algebras (Aij), each with a unit lij and a state (pij, and let 
(Tijj, vTjj, (^jj) be the associated GNS triples, so that (pij{a) = (ji,jia)ii,j,^i,j} for any 
a e Aij. For any let Xij be the *-representation of Aij given by 

Xi,j{a) = Vij{7iij{a) ® In{i,j))V*j for a e Alj, 

where In(i,j) denotes the identity on H{i,j). Note that these representations are, in 
general, non-unital. In fact, 

where and Sjj are canonical projections in B{7i) given by 

ri,j = P-Hiij) and Sij = PK^ij) 

where JC{i,j) = T-ti j®H{i,j). For given the projections j and Sj j are orthogonal 
and their sum is the canonical projection onto the subspace of H onto which \{Aij) 
acts non-trivially. 

The Xij remind the representations Aj of free probability, but the corresponding 
operators Ajj(a) have larger kernels. Using Ajj's, we shall define product representations 
on 

:= UijAij, 

the free product without identification of imits, equipped with the unit 1^, and products 
of states which are analogs of the free product representation and the free product of 
states, respectively. 

Definition 2.3. The matricially free product representation ttm = is the unital 

*-homomorphism A : ^ ^ B{T-L) given by the linear extension of 

A(l^) = 1 and A(aia2 . . . a„) = Xi-^j^{ai)Xi.^j.^{a2) . . . Xi^j^{an) 

for any ak e Ai^^j^, k = 1, . . . ,n, with (hji) ^ {12, j2) ^ ■ ■ ■ ^ (injn)- The associated 
state (/? = : ^ ^ C is given by 

^{a) = <7rM(a)^,0 
and will be called the matricially free product of (ipij). 

Basic properties of the product state <fi are collected in the propositions given below. 
Roughly speaking, they show that this state (on the free product of C*-algebras with- 
out identification of units) has similar properties as the free product of states (on the 
free product of C*-algebras with identification of units) except that the units of these 
algebras act as units only on 'matricial' tensor products and otherwise they act as null 
projections. 

For that purpose, it will be useful to introduce the following sets of indices: 

An = {((^1,^2), («2, is), {in, in)) ■ ^2) ^ («2, is) ^ ■ ■ ■ ^ {in, in)} 

and A = U^_iA„. Clearly, the set A encodes the matricial, freeness and diagonal 
subordination properties. Finally, X stands for the unital subalgebra of A generated by 
the units lij. 
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Proposition 2.2. Let ip be the matricially free product of states (fij) and let e 
Ai^,j^, where k e [n] and (ii, ji) ^ . . . ^ {i 

nt jn) ■ 

(1) // ttfc G Ker ipi^^j^ for k e [n], then f{aia2 ...a„) = 0. 

(2) // ar = l^j, and e Ker ipi^j^ for r < m ^n, then 

ip{ai...ar-iar+i...an) if {{ir, jr), ■ ■ ■ , {^n, jn)) ^ ^ 

otherwise 



ip{ai ...an) 



(3) For any a e A, Ui, U2^T and i, j e I , it holds that 

(p{uiau2) = ip{ui)(p{a)if{u2) and (fiiUj) = Sij. 

(4) The restriction of to Ajj is (pjj for any j e I. 

(5) The mixed moments (p{aia2 ■ . . an) are uniquely expressed in terms of mixed 
moments of products ofak's in the states fi/.j^.. 

Proof. If ak e Ker(/9j^jj. for k e [n], where 7^ ... 7^ {in,jn), then it follows 

from the definition of the Xij that 

(1) nM{ai...an)C = if {{ii, ji), . . . , {in, jn)) ^ A 

(2) nMiai...an)ienlj^®...®Hlj^ if ji), . . . , g A. 

In both cases we obtain a vector orthogonal to ^ on the RHS, which proves (1). Sup- 
pose now that the assumptions of (2) hold. If {{ir,jr), ■ ■ ■ , {in,jn)) ^ A, then ^irjriUrjr) 
acts as a unit onT-G , ® . . . ® Ti^ ■ by the definition of the representations Aj 
On the other hand, Xirjri'^irjr) kills any simple tensor beginning with h e 'H^^^^^^^^ if 
jr ir+i or {{ir+ii jr+i)i ■ ■ ■ , {inijnj) ^ A sincB Vi^j^ docs, which completes the proof 
of (2). In turn, (3) follows from the action of the A(lij) onto ^. That agrees with 
(Pjj on Ajj for any j e I follows from the action of the Xjj{a), a e Ajj, onto ^, namely 
7rM(fl)C = + '/'jjC^)^! which gives (4). Finally, (5) is a consequence of (l)-(2). 



In a similar way we can define states associated with other unit vectors from Ti,. For 
our purposes, we will need states associated with unit vectors Cj e Ti.'j j which are in 
the ranges of '^j,j{Ajj), where j e /, respectively, namely (pj : A ^ C defined by the 
formulas 

ipj{a) = {nM{a)ej,ej). 

These states will be called conditions associated with </? and will be used for computing 
traces. They have similar properties as (p as the proposition given below shows in more 
detail. The main difference is that (pj satisfies the condition of freeness type only for 
indices which satisfy {j,j) 7^ {ii,ji) 7^ ... 7^ {imjn) ^ (JjJ)- Moreover, ipj\I is quite 
different than ip\I due to different normalization conditions. 

Proposition 2.3. Let (pj, where j e I, be the conditions associated with ip = 

and let let ak e A^j-fc, where k e [n] and {j,j) 7^ {ii,ji) 7^ ... 7^ (injn) ^ ij,j)- Then 

(1) If ak G Ker for k e [n], then pj{aia2 . . . an) = for each j. 

(2) // ttr = Iv J, and a^ e Ker Pi^j^ for r < m ^ n, then 

( . _ ( (Pj{ai. . .ar-lttr+l. . .an) if {{ir,jr),---,(in,jn))^^ 

iPj[a^...an) - < Q otherwise 
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(3) 



For any a e A, Ui,U2 and k e I, it holds that 

(Pj{uiau2) = ipj{ui)(pj{a)ipj{u2) and ^j{li,k) = ^j,fc- 



(4) The restriction of ^pj to Aij is (fij for i ^ j. 

(5) The mixed moments (pj(aia2 ■ ■ ■ an) are uniquely expressed in terms of mixed 
moments of products ofak's in the states ^Pi^.,jf.- 

Proof. The proof is similar to that of Proposition 2.2. The main difference concerns 
(3) and (4). In this context, notice that the unit vectors ej play the same role with 
respect to the action of the Xi,j{a) for any i ^ j as plays with respect to the action 
of 7rij(a), where a e Aij, and thus (pj{a) = (Xij{a)e-j,ej) = ('7^i,j{a)Ci,j,Ci,j} = ^iM, 
which gives (4). ■ 

Remark 2.1. Note that states f and (</7j) share together the property of extending 
the array of states (<^ij). Thus, ip extends the diagonal states Lpjj for all j, but it does 
not extend the non-diagonal states for i j since TTM{o)i = for any a 6 A;/,. In 
turn, Lpj extends ipi^j for i ^ j, but it does not extend ipjj or any of the other diagonal 
states. This is a natural consequence of differences in the definitions of the diagonal 
and non-diagonal partial isometrics. Moreover, it is clear that for each ipj there exists 
hj G Aj^j n K.erip such that 



for any w e UijAij. Therefore we can reduce computations of mixed moments in the 
state (fj to computations of mixed moments in the state (p. 

Finally, let us denote by the unital commutative *-subalgebra of S(7i) generated 
by the A(lij), where i,j e I. By abuse of notation, will also be denoted by l^j. 

Observe that these projections are not mutually orthogonal and that is why it is often 
convenient to use their subprojections r^j and Sjj. 



Of special importance is the subspace of the matricially free Fock space, called the 
'strongly matricially free Fock space', in which the diagonal Hilbert spaces appear only 
at the end of tensor products. The main reason is that it is related to both free and 
monotone Fock spaces. We also study the associated product states which can be viewed 
as direct generalizations of both free and monotone products of states. 

Definition 3.1. By the strongly matricially free Fock space over Ti. := {Ti-ij) we under- 
stand the subspace of A^(7Y) of the form 



3. Strongly matricially free products 




m=l ii^.-.^im 
,...,n7TiGN 



with the canonical inner product. 



A justification for the word 'strong' is that in this case the words iiZ2 ■ ■ - im which 
label the tensor products in the above definition satisfy ii 7^ i2 7^ ■ ■ ■ 7^ im- 
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Example 3.1. The simplest space of this type is associated with a 2-dimensional square 
array Ti.. Then 

00 

n{H) = 07^('")(^)> 

where the first few summands are of the form 

n^'Hn) = 7^1,1® 7^2,2 

n^^\n) = w^?®wg®(^i,2(8)W2,2)®(H2,i(8)7ii,i) 

7^(3) (7^) = 7^®3 ^ ^ (^^^^ ^ ^®2^ ^ ^^^ ^ ^ ^®2^ ^ (^®2 jg, ^^^^^ 

®(7^g ® 7^2,2) ® {'Hl,2 ® ^2,1 ® ® (7^2,1 ® :^1,2 ® ^2,2), 

etc. In contrast to M.{i-C), we do not have tensor products like 7-^2,2 ®'H2,i ®'Hi,i and 
7^1,1 ® 7^1,2 ® 7-^2,2 in the summand of the third order. 

Remark 3.1. For a given array of Hilbert spaces 7Y = (7ijj), we have inclusions 

n{n) ^ M{n) ^ n®H,,i) 

which, in most cases, are proper. Moreover, if we have a square array and Tiij = Tii 
for any G /, where {'Hi)iei is a family of Hilbert spaces, then there is a natural 
isomorphism 

iel 

Since 7^®% 07{^;^ ® . . = ® ® ■ ■ • ® ^t™ ^r any Zi, Z2, . . . , g /, 

ni, . . . Urn, m G N. Similarly, if we have a lower-triangular array and 7iij = Tii for any 
i ^ j, then IZiTi) is isomorphic to the monotone Fock space. 

Moreover, as expected, there is a product of Hilbert spaces related to the strongly 
matricially free Fock space, and an analog of Proposition 2.1 holds. 

Definition 3.2. By the strongly matricially free product of {Ti-ijj^ij) we understand 
the pair (Q,^), where Q is the subspace of 7i of the form 

00 

We denote it {0,0 = 

If we consider a family of unital C*-algebras {Ai)iei, each equipped with a family 
of states {ipij)jej, then we can look at this product space as follows. If {Hi,j,T^i,j,ii,j) 
is the GNS triple associated with the pair {Ai,tfiij), then the Hilbert space Tijj (as 
well as the corresponding state and representation) taken as the representation space 
for the algebra Ai at some given tensor site depends on the algebra Aj represented at 
the following tensor site on the space 'Hj^k for some k. It is worth noting that in this 
framework we can assume that for fixed i g 7 all vectors g /, are identified since 
we can take the tensor product of Hilbert spaces ®j'Hi^j and set = for each 
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i e I. It is not hard to see that in this framework our model is related to freeness with 
infinitely many states [8,9]. 

The construction of the product state is similar to that in the usual case. The only 
difference in all definitions is that the sets A„ are replaced by 

r„ = {((^1,^2), (^2, ^3), ■ ■ ■ , (in, in)) : il 7^ i2 ^ ■ ■ ■ ^ in] 

and their union A by F = (Jl^^^ r„. Note that the conditions which define the r„ are 
stronger than those which define the A„ and therefore all objects constructed in the 
strong case are obtained from the standard ones by a projection-type operation. 

The partial isometrics in the strong case, denoted by Wij, remind the Vij except 
that they refer to Q rather than H. In particular, the diagonal partial isometrics 
'■ '^3,3 ® ^UJ) Q are given by 

and nl®^^nl 

where S{j,j) = for any j, whereas the non-diagonal partial isometrics Wjj and the 
associated subspaces G{i,j) are similar to the T-l{i,j). 

Definition 3.3. Let Pij be the *-representation of Aij defined as in the standard 
case, with Vij replaced by Wij for any The corresponding strongly matricially free 
product representation TiR = and strongly matricially free product of states 

are defined in terms of the pij as in the standard case. 

The basic results on the strongly matricially free product of states are similar to those 
in Propositions 2.2-2.4 and therefore we state them in an abbreviated form without a 
proof. 

Proposition 3.1. Let Lp be the strongly matricially free product of states {fij) and 
let ipj be the state associated with any unit vector from H'j j which is in the range of 
Pj,j{Ajj), j G /. Then the statements of Propositions 2.2-2.3 remain true, with A 
replaced by V. 

As we mentioned earlier, one of the advantages of using the strong structures is that 
they are straightforward generalizations of those in free probability (in the case of square 
arrays) and monotone probability (in the case of triangular arrays) . 

Remark 3.2. Using the strongly matricially free product of states, we can obtain the 
free product of states if we take a square array. Namely, if Ai^j = Ai and (pij = (pi for 
all i,j G /, then, on the level of Hilbert spaces, we have Tiij = Hi and = for all 
i,j G /. Thus 

(g,e) = *.e/(^.,6) 

and pi{a) := '^jPijia), a e Ai, understood as the strong limit, extends to a (unital) *- 
representation of Ai on Ti.. Then p = Ui^iPi agrees with the free product representation 
of UieiAi on Q and therefore the corresponding strongly matricially free product of 
states gives the free product of states. If (pi^i = ipi and ipij = ipi for any i 7^ j, we obtain 
the conditionally free product of states. 

Remark 3.3. The monotone product of states [17] is obtained from the strongly ma- 
tricially free product of states when we take a triangular array T, with Aij = Ai and 
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= for all {i, j) e T I X I (it can be obtained from the square array by setting 
Aij = for i < j), where / is a totally ordered set. Then Hij = Hi and = for 
i ^ j and 

the expression on the right-hand side being the monotone product of Hilbert spaces. 
Moreover, rj(a) := Yjjs^iPhA'^) extends to a (non-unital) *-representation of Ai on 
such that r = Ui^jTi agrees with the monotone product representation of Ui^iAi on 
Q and therefore the corresponding strongly matricially free product of states gives the 
monotone product of states. 

4. Matricial freeness 

Guided by the notion of the matricially free product of states, we shall introduce now 
the associated concept of independence called 'matricial freeness', and closely related 
to it, 'strong matricial freeness'. They involve arrays of noncommutative probability 
spaces and to some extent they remind models with many states [6,8,9], but they cannot 
be reduced in a natural way to any of these (freeness with infinitely many states has 
some non-empty intersection with 'strong matricial freeness' and conditional freeness 
is its special case). Moreover, we will study discrete matricially free Fock spaces, both 
standard and strong. 

Let ^ be a unital algebra with an array (Aij) of not necessarily unital subalgebras 
of A. Assume that each Aij has an internal unit Ijj and that the subalgebra X of ^ 
generated by all internal units is commutative. Further, let (f he a distinguished state 
on A such that = Sij for any i,j and let {tpj : i e /} be a family of additional 

states on A such that (pj(li^k) = ^j,k for ciny i, j, /c, which will be called conditions. Here, 
by a state on A we understand a normalized linear functional (if .4 is a *-algebra, we 
require this functional to be positive). 

In this situation it is convenient to form an array (<^ij) of states on A by 

^Pj^j = (/? and ^i^j = ^pj for i ^ j 
which will be said to be defined by the state (/? and the conditions ipj. 

Definition 4.1. Let (v^jj) be defined by the state tp and the conditions pj. We say 
that (Ijj) is a matricially free array of units associated with {Ai,j) and {pi^j) if 

(1) Lp{uiau2) = f{ui)p{a)ip{u2) for any ae A and ui,'U2 e X, 

(2) if Ofc G ^fcjfc n Ker(^jj.jj^,, where r < k ^ n and r < n, then 

, , j ip{aar+i...an) if ((v, >), . . . , (4, Jn)) e A 

c^(al,,,,-,a.+i...a„)-| ^ otherwise 

where a e ^ is arbitrary and {ir,jr) ^ • • • ^ (^m jn)- 
The array (Ijj) is called a strongly matricially free array of units if A is replaced by F. 

The above definition enables us to define the concepts of matricial freeness and its 
strong version called strong matricial freeness. They both bear some resemblance to 
freeness, but the main difference is that the identified unit in the context of freeness 
is replaced by the (strongly) matricially free array of units. In fact, we will see later 
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that it is the strong matricial freeness which can be viewed as a direct generalization 
of freeness. 

Definition 4.2. We say that (Aij) is matricially free with respect to (i^ij) if 

(1) for any e Keiipi^j,^ n Ai^j^, where k e [n] and {ii, ji) 7^ . . . 7^ {in,jn), 

(p{aia2 . . . a„) = 

(2) (Ijj) is a matricially free array of units associated with (Aij) and {'^ij). 
In an analogous manner we define strongly matricially free arrays of subalgebras. 

Definition 4.3. The array of variables (ajj) in a unital algebra A will be called 
{strongly) matricially free with respect to the array (^ij) defined by a state (p and 
the conditions ipj if there exists an array of (strongly) matricially free array of units 
(Ijj) in A such that the array of algebras (C[ajj, Ijj]) is (strongly) matricially free with 
respect to {^ij). 

Using the above definitions, we can compute mixed moments of arbitrary matricially 
free random variables. Namely, writing = a° + fik,jki^k)hk,jk each k e [n], we 
obtain a recursion 

ip{ai...an) = 2 (^i,j,(afe)(^(a;...li,j, ...<) 

IsCfcsCra 

+ Ya '^^fe.ife i'^k)<^ii,3i iai)M ■ ■ ■ hkJk ■ ■ ■ 

l^k<l^n 

+ 

+ ^njl (Ol) ■ ■ ■ ^in,jnM<f{hl,h ■ ■ ■ K,jn) 

for any (^l,jl) 7^ ... 7^ {in,jn)- It is easy to see that that the mixed moments on 
the right-hand side, written here in a slightly simplified manner, reduce to moments of 
orders smaller than n. Note that these depend on (fi\I in an essential way. That is why 
the states tp and the ipj of Section 2 give different mixed moments, although they have 
similar properties. 

Remark 4.1. If ^ is a unital *-algebra, then, in addition, we require that the func- 
tionals ipij are positive, the Aij are *-subalgebras and the lij are projections. Then 
an array of variables (a^j) will be called *-{strongly) matricially free if the array of 
*-algebras (C(aij, a*j, Ijj)) is (strongly) matricially free. 

Example 4.1. Let us assume that we have a two-dimensional square array of variables, 
namely ai,i = a, ai,2 = a',a2,i = b',a2,2 = b, which are strongly matricially free with 
respect to the array {(fij) defined by ip and the pair (ipi, (P2). Using the above recursion 
and then Definitions 4.1-4.2, we obtain 

ifiiab'ab) = ^(6)^i(6')(<^(a')-<^'(a)), 
(f{abab) = ip^ {a)(p^ {b) , 
<p{aba'b) = ip{a)ip2{a'){<p{b'') - cp^b)). 

Note that if we assume that 9^1 = </?2 = i^, then the sum of these moments gives (p{abab) 
for a, b conditionally independent with respect to {(p, ip). In fact, the above sum is equal 
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to ip{{a + a'){b + b'){a + a'){b + b')). In turn, if we take the lower-triangular subarray 
and take ^pi and (/^-distributions to be equal, the sum of the first two moments (the 
third one is zero since there is no a') gives Lp{abab) for a, b monotone independent with 
respect to (f Similar agreements hold for mixed moments of higher orders, which is in 
agreement with Remarks 3.2-3.3. 

Definition 4.4. By a discrete matricially free Fock space we understand M. = M.{7i), 
where Tii^j = Ccjj for any (i, j) G N, and the array (cjj) forms an orthonormal basis of 
some Hilbert space. In an analogous manner we define the discrete strongly matricially 
free Fock space TZ = TZ{TC). 

Both Ai and TZ are subspaces of the discrete free Fock space ^ Ccij) and thus 

allow for the canonical action of free creation and annihilation operators. In order to 
specify the arrays of units, let us distinguish two types of their subspaces. 

In the case of Ai these are 

(1) A4(i,j), spanned be simple tensors which begin with Cj^k for some k, where 
(i) ^ and, in addition, by f2 if i = j, 

(2) IC{i,j), spanned by vectors which begin with e^j, where i,j are arbitrary, 

and the direct sum IC{i,j) is the subspace of Ai onto which the *-algebra 

generated by i{eij) restricted to M. acts non-trivially, where i{eij) denotes the canonical 
free creation operator associated with vector e^j. 

The canonical projections onto such direct sums are natural candidates for the ma- 
tricially free units Ijj and therefore we set 

U,j ■= PK.{i,j)®M{i,j) and = i{eij)lij 

for any The adjoint of £ij will be denoted by Note that the projection Ijj is an 
internal unit in the *-algebra Aij = C{£ij,£*j) and £*j£i,j = for any However, 
in the algebra C{iij,i*j : i,j e N) there are more relations as the proposition given 
below demonstrates. 

Proposition 4.1. In the algebra C(£ij,i*j i,j ^ N) the following relations hold: 

(1) ^Ijkj = hj for any ij, 

(2) £lj£k,i = 0, iijik.i = and li,j£k,i = whenever {i,j) 7^ {k,l) and j 7^ k, 

(3) li,/j,k = ^j,k for any j, k. 

Proof. We omit the elementary proof. ■ 

In the case of TZ, we proceed in a completely analogous fashion and distinguish the 
following subspaces: 

n{i,j)=M{iJ)nTZ and C{i, j) = K{i, ]) nil 

which lead to the definitions of the strongly matricially free array of units and of the 
creation operators, respectively, 

= Pc{i,j)®n{i,j) and kij = i{eij)lij 
where, slightly abusing notation, we use the same symbols for the units as before. 
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Finally, we specify the states: will denote the vacuum state associated with fi, the 
</7j will be the states associated with the e^j, j e N, and ((^ij) will be the array defined 
by if and the ipj. Here, the same notation is used for states on B{Ai) and B(TZ). 

Proposition 4.2. The array of *-subalgebras Aij = C(£,^,, ^f -) of B{M), where i,j e 
N, is matricially free with respect to {'^ij)- The array of *-subalgebras Bij = C(kij, k*j) 
of B(TZ), where i,j e N, is strongly matricially free with respect to {(fiij)- 

Proof. The proof is similar to that of Voiculescu for the discrete free Fock space given 
in [21]. We shall look at the case of (Aij) since the case of (Bij) is analogous. Each 
algebra Aij is spanned by operators of the form 

^l/ij, where p + q > 0, 
and the projection Ijj. However, the corresponding moments vanish: 

for any i,j since p + q > 0. Moreover, j) = 1 for any Therefore, in order to 

show that condition (1) of Definition 4.2 holds, it is enough to show that 

whenever (ii, ji) ^ . . . ^ {in,jn) and pi + qi > 0, ...,;?„ + g„ > 0. The same argument 
as in [21] allows us to reduce the proof to the case when (?i = . . . = g„ = 0, which 
implies that pi > 0, . . . ,Pn > 0. But then the moment clearly vanishes. This proves 
condition (1) of Definition 4.2. Condition (2) follows easily from the definition of the 
projections l^j in view of the relations given in Proposition 4.1. This completes the 
proof. ■ 

Example 4.2. For i,j e I, let Gij = F{1) be the free group on one generator gij 
with unit e^j, and let \ij be the corresponding unitary operator on the space P{Gij) 
given by Xij{g)5{h) = S{gh), where {S{g) : g e Gij} is the canonical basis. Consider the 
subspace 1^ of l'^{*i,jGij) spanned by vectors of the form S{g), where g is either the unit 
e of the free product *i^jGij, or a product of the form gig2 ■ ■ ■ gm, where gk e Ci^.i^+i • = 
^ikik+i\{^ik,ik+i} foi' 6ach k, with (^l,^2) ■ ■ ■ («m, Wi) and i^+i = im- The space 
l\f can be viewed as the space of square integrable functions on the (non-existent) 
'matricially free product of groups'. Let Ij j denote the projection from l^{*ijGij) onto 
the subspace of spanned by vectors S{g), where g begins with an element from G^j 
or for some k or, in the case of i = j, also g = e. Define 

for g e Gij. Then the array (Aij), where Aij is the *-subalgebra of B{l'^{*ijGij)) 
generated by K,j{9i,j) and lij, with the standard involution, is matricially free with 
respect to the array where the diagonal states coincide with (p{.) = (.(5(e), (5(e)) 

and the non-diagonal ones in the j-th column coincide with <^j{ ) = {■S{gjj),S{gjj)}. 

Example 4.3. In the above example take an n-dimensional square array of copies of 
F(l) and denote their generators by gij, where i,j e [n]. For each natural n we form 
a tree (a subtree of the homogenous tree ]Hl2„2) which corresponds to the 'matricially 
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free product of free groups'. Suppose the root e corresponds to the 'father'. We 
distinguish 'sons' and 'daughters' in each 'generation' which correspond to the left 
action of Qj^j or gj^, and gi^j or g^ j, respectively, where i ^ j. The rules of drawing the 
tree follow from matricial freeness and are the following: each 'son' has 1 'son' and 2n—2 
'daughters', whereas each 'daughter' has 2 'sons' and 2n — 1 'daughters'. Therefore, 
'sons' and 'daughters' correspond to vertices of valencies 2n and 2n + 2, respectively. In 
Figure 1 we draw such a tree for n = 2 (black and empty circles are assigned to 'sons' 
and 'daughters', respectively). If we make an additional assumption, for instance that 
'daughters' cannot have 'sons' (this fact corresponds to strong matricial freeness, where 
diagonal generators kill words beginning with the non-diagonal ones), we recover IH[2„ 
of free probability (cf. Remark 3.1). 

Example 4.4. Let TZ{H) be the discrete strongly matricially free Fock space and let 
TZ{J-L) = ^(©jCe,) be the natural isomorphism of Remark 3.1, where {cj : j e N} is 
an orthonormal basis of some Hilbert space. If {wi^j) is an infinite matrix with non- 
negative parameters p and q above and below the main diagonal, respectively, and I's 
on the diagonal, then this it is easy to see that the (p, g)-creation operators studied in 
[15] can be identified (with the use of this isomorphism) with the strongly convergent 
sums 

j 

where i e N and the (p, g)-annihilation operators are their adjoints. A similar approach 
can be applied to square arrays of arbitrary Hilbert spaces and *-representations, which 
leads to some notion of '(p, g)-independence'. Moreover, it can be carried out for more 
general matrices (wij) within the framework of the strong matricial freeness, which 
generalizes notions of independence of this type. 
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5. Traces 

In this Section we introduce some real-valued functions on the set of non-crossing pair 
partitions. These functions are obtained by computing traces of a square real-valued 
matrix. We will assume later that this matrix has non-negative entries which represent 
variances of probability measures on the real line M.r and we will demonstrate that the 
functions introduced in this Section describe the asymptotics of matricially free random 
variables in central limit theorems. 

Let AfCm denote the set of non-crossing partitions of the set [m], i.e. if tt = 
{7ri,7r2, . . . jTTfe} G HCm, then there are no numbers i < p < j < q such that i,j e tt^ 
and p,q e Tig for r ^ s. The block tt^ is inner with respect to tt^ if p < z < g for any 
i e TTr and p,q e Tig (then vr^ is outer with respect to tt^). It is clear that if tt^ has outer 
blocks, then there exists a unique block among them, say tt^, which is nearest to tt^, i.e. 
if another block, say tt^, is outer with respect to tt^, then we must have a < p < 6 for 
any a,b e TTt and p e tt^. Then the pair (tt^, tt^) is called the nearest inner-outer pair of 
blocks. 

Let MCCm denote the set of non-crossing covered pair partitions of [m], by which we 
understand the subset of MCm consisting of those partitions in which 1 and m belong 
to the same block (if m = 1, we understand that the partition consists of one block). 
In terms of diagrams, all blocks of tt e HCCm, where m > 1, are covered by the block 
containing 1 and ni. We denote by A/'C^ and AfCCl^ the sets of non-crossing pair 
partitions of [m] and non-crossing covered pair-partitions of [m], respectively, and we 
set MC' = [jZ^,^fCl and AfCC' = [jZ^^MCCl. 

It is easy to see that each tt g N'Cm can be decomposed as 

(5.1) TT = TT^^) U 7r(2) U . . . U TT^P'^ 

where 7i^^\ .... tt^^^ are non-crossing covered partitions of subintervals h, I2, ■ ■ ■ , Ip of 
[m] whose union gives [m]. By a partition of a set / consisting of r elements we 
understand the corresponding partition of [r]. 

On the other hand, each tt e AfCCm can be decomposed as 

(5.2) TT = 7r(°) u 7r(i) u . . . u tt^'^ 

where tt^'^^ is the block containing 1 and m and TT^^\7r^'^\ . . . , tt^'^^ are non-crossing cov- 
ered partitions of subintervals /i, I2, . . . ,Ir of the set {2, . . . , m — 1}. 

Consider now a square real-valued matrix V = (vij) e M„(R), where n e N u {co}. 
The usual trace and the normalized trace will be denoted 

n 1 ™ 

Tv{V) = Y,Vj,j and tr(y) = - 

3=1 j=i 

respectively. For finite n, we define the 'diagonalization mapping' 

T : M,(M) ^ D^M), t{V) = diag(2^,,i, • • • ,X!^i,n) 

j j 

where -Dn(l^) is the set of square diagonal real- valued matrices of dimension n (by abuse 
of notation, the same symbol r is used for all n). In other words, r computes the sum 
of all elements of V in each column separately and puts this value on the diagonal. 
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Figure 2. Colored partitions with conditional blocks. 



Using the above trace operations, we shall define two real-valued functions on the 
set denoted v and vq, associated with given V e M„(R). Although there is a 

close similarity between these functions when restricted to AfCC^ (in particular, they 
are defined as traces of certain matrix- valued quasi- mulitplicative functions), note that 
they are extended to MC^ in two different ways. 

Definition 5.1. For a given matrix V e M„(R), we define a mapping from jVC^ to 
I?„(]R) by assigning to each tt g AAC^ the matrix V{Tr) by the following recursion: 

(1) if 71 consists of one block, we set V{7r) = t{V). 

(2) if TT G J\fCC^ consists of more than one block, then 

(5.3) V{7r) = r{V{7r^^^) . . . ^(Tr^'-))^) 

according to the decomposition (5.2), 

(3) if 7^GA/■C^ then 

(5.4) V{7t) = y (7r(^))y (tt^")) . . . V{tt^p^) 

according to the decomposition (5.1). 
Let t) : A/'C^ ^ M be the function defined by v{n) = tr(V(7r)). 

Example 5.1. For some V e M„(]R), consider three partitions given in Fig. 2. Color 
each block by a number from the set [n]. Computation of the corresponding values of 
the function v gives 

v{n) = tv{riT{V)V)) = ^J]v,,vj,, 

virj) = tT{riT{V)V)TiV)) = ^Y^v,,vj,iVk,i 

i,jyk,l 

viC) = tv{riT{V)riV)V)) = ^ v,,kVj,kVk,i 

i,j,k,l 

One can see that to each nearest inner-outer pair of blocks (vTr, tt^) we assign the matrix 
element Vp^q, where p and q are the colors of Hr and tTs respectively. Moreover, if a 
block does not have any outer blocks and is colored by q, then we assign to it the 
matrix element Vg^t, where t is assumed to be the same for all such blocks (one can 
imagine that we have an additional 'conditional block' colored by t which covers all 
other blocks). At the end we sum over all colorings. 
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The function Vq is defined on the set AfCC^ in a very similar manner, except that we 
replace the right multiplier V in (5.3) by its main diagonal 

Vo := diag(^;i,i, . . . ,^;„,„) 

which corresponds to changing only the contribution of the covering block. Then we 
extend vq to all of jVC^ by multiplicativity of vq. 

Definition 5.2. For a given matrix V e M„(R), we define a mapping from AfCC'^ to 
D„(R) by assigning to each tt g J\fCC^ the matrix Vo(7r) by the following recursion: 

(1) if TT consists of one block, we set Vo(7r) = Vq 

(2) if TT G MCC^ consists of more than one block, then 

(5.5) yo(7r) =T(y(7r(i))...y(7r('-))yo) 

according to the decomposition (5.2). 

Let : MC^ ^ M be the function defined by vq{t\) = Tr(V'o(7r)) for tt g J\fCC^, extended 
to A/'C^ by mulitplicativity Vo(n) = i'o(7r(^)) . . .Vo{n^^^) according to the decomposition 
(5.1). 

Example 5.2. Let us compute the values of vq corresponding to the partitions in Fig. 2. 
We have 

Vo{n) = Tr{T{V)Vo)=Y,%jVj,j 

Mv) = Tr{T{V)Vo)Tr{Vo) = ^ Vi,jVjjVk,k 

Vq{,C) = Tr(r(r(y)r(y)Vo)) = 2 ^^i,fe^i,fcVfc,fe 

Note that the main difference (apart from normalization) between vq{t:) and f (tt) con- 
cerns the blocks which do not have outer blocks. Here, if such a block is colored by g, 
we assign to it the matrix element Vq^q and we do not use 'conditional blocks'. 

Below we shall prove a lemma which gives explicit combinatorial formulas for v{tt) 
and vo{ti). If tt = {tti, . . . ,7rfc} g MC\^,, we color the blocks of vr and the conditional 
block TTo by the set \n\. The coloring of the partition tt u {ttq} can then be identified 
with the mapping / : [A;] u {0} [n], where f{i) is interpreted as the color of tTj. The 
pair (tt, /) can then be interpreted as the colored partition associated with tt and the 
mapping /. The set of all colorings of tt by the set [n] will be denoted Fn{'K). 

For a given matrix V = (vij) e Mn(M) and given colored partition (tt, /), where tt is 
a non-crossing pair partition consisting of k blocks and / is the coloring of tt, we can 
now define (slightly abusing notation) a natural mapping 

i;o:{(7ri, /),..., (7rfe,/)}^M, 

by the following rule: 

voiT^pJ) = Vi,j if f{p) = i and f{a{p)) =j 

where a{p) is the index of the nearest outer block of TTp if tt^ has an outer block, and 
otherwise a{p) = p. In other words, to each block we assign the matrix element whose 
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first index is the color of that block, and the second index - the color of its nearest 
outer block. 

The multiplicative extension of this function to the class of non-crossing colored 
partitions gives 

Vo{t^, f) = VoilTi, /)Vo(7r2, /) . . . Vo{Trk, /). 

In a similar way we define functions 

for any j e [n], where Vj{7rp,f) = Vij if f{p) = i and f{(Tj{p)) = j where (Tj{p) is the 
index of the nearest outer block of tt^ if tt^ has an outer block, and otherwise we put 
(7j(p) = j. In the latter case, j is interpreted as the color of the conditional block. 

Lemma 5.1. For any V e M„(R), where n g N, and tt g MC'^ , it holds that 
^o(7r) = 2 Vo{'kJ) and v{'k) = - ^Uj(7r,/) 

where the summation over j corresponds to all colorings of the conditional block. 

Proof. We provide an induction proof for the function v (the proof for ?,'o is similar) . 
The main induction step will be carried out on the level of the (diagonal) matrices 
V{7r). We claim that its diagonal entries are of the form 

for any tt g MCC'^ and qe \n\. In view of (5.4), the required formula for v{ti) is then a 
straightforward consequence of the claim. Of course, if tt consists of one block, then 

and thus our assertion easily follows. Assimie now that tt has k > 2 blocks and suppose 
the assertion holds for non-crossing covered partitions which have less than k blocks. 
Since tt has a decomposition of type (5.2), the assertion holds for tt^^^tt^^^ . . . ,7r^^) in 
this decomposition. We know that the matrix assigned to tt has the form (5.3) and 
therefore the product of (diagonal) matrices corresponding to these subpartitions has 
(diagonal) matrix elements of the form 

{V{'k)U = {t{WV)\,, = Y,W,,v,,, 

3 

where q is the color of the conditional block of tt and 

Wi,j= 2 V,{'K^'\f{) ... 2 V,{'K^^\fr) 
/ieF„(7r(i)) /reF„(7r('-)) 

by the inductive assumption. Now, since the blocks oi tt^^\ . . . ,tt^'^^ are colored inde- 
pendently, we have 
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for a uniquely determined coloring / of the blocks of tt in which j can be interpreted as 
the color of tt^^^ since vr^") covers vr^^^ . . . ,7i^^'> and q can be viewed as the color of the 
conditional block of tt. This and the above formula for Wjj gives the desired formula 

feFnin) 

which proves our claim and thus completes the proof of the theorem. ■ 

Under suitable assumptions on V, there is a simple connection between functions v 
and Vq if Vq = ain/n, where is a unit matrix and a is a positive number. We shall 
express this relation in terms of the corresponding formal Laurent series, which turn 
out to be Cauchy transforms of (compactly supported) probability measures on the real 
line associated with appropriately constructed random variables, 

Proposition 5.1. Let V e M„(R) be such that Vq = aIn/n, where a > 0. If G{z) = 
XifeLo '^2fe-2^^'^^^ and Gq{z) = XlfeLo ^sfc-^^^'^^^ are formal Laurent series, where 

a2k = Yj '"^^^ ^2fe = X; ^o(^) 

and both v{7r) and vo{Tr) are associated with V , then Gq{z) = Xjiz — aG{z)). 
Proof. Observe that in the case when Vjj = a/n for all j e [n], we have 

r N ^o(7r') 

vItt) = — — - 
a 

for any n e A/'C^, where the partition n' e J\fCC^j^2 i'' obtained from tt by adding to tt 
the block that covers all blocks of tt, say {0,m + 1}. This leads to 



2m 



a "—t, az^ 



where C{z) = X/m=o C2m-2^^"^ and C2m = YjweMcc'^ ^o(7r). Xow, using the multiplicativity 
of vo, we obtain 

00 00 -. 

B{z) := Y h,mz^- = 1 + X {C{z) - \T = 



m=0 m=l 



which leads to 



where 



M^) = ^ ( 1 - -Jr?\ and Go{z) ^ 



az^ V B{z) J ' z-aG{z) 



G(z) = -A (-] and Go(z) = -B (- 
z \z I z \z 



which completes the proof. 
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6. Random pseudomatrices 

In this Section we will study the asymptotic behavior of random pseudomatrices for 
two sequences of states: the sequence of distinguished states, with respect to which our 
variables will be matricially free, and the sequence of traces which are normalized sums 
of conditions in the definition of matricial freeness. Under suitable assumptions, we will 
later obtain two types of central limit theorems: the standard central limit theorem as 
well as the tracial central limit theorem for matricially free random variables, the latter 
being related to random matrix models. 

Let {An)nen be a sequence of unital *-algebras. For each n, let {Xij{n))i<^ij<^n be an 
array of self-adjoint random variables in An and let {(j)ij{n)) be an array of associated 
states on An- We will say that the moments of the Xij{n) are uniformly bounded with 
respect to [(f)i^j{n)) if, for all natural m, there exists ^ such that 

(6-1) \^iAn){X-{n))\ ^ 

for all n G N and i, j e [n]. 

Assume further that for each natural n the array (Xj ,, (n)) is matricially free with 
respect to the array {4>ij{n)) defined by a distinguished state 0„ and the conditions 
ipn,j, where 1 ^ j ^ n. We are going to study the asymptotic behavior of random 
pseudomatrices 

n 

(6.2) S{n) = 

with respect to two types of states: 

(1) distinguished states 0„, 

(2) convex linear combinations of conditions 

1 

(6.3) ■= - I] ^".i 

which play the role of traces in the context of random pseudomatrices. 

This will lead to two types of central limit theorems for matricially free random 
variables: the 'standard CLT' reminding the CLT for free random variables [20] and 
the 'tracial CLT' reminding the limit theorem for random matrices [22]. Thanks to the 
'matricial property' and the 'diagonal subordination property' of the matricially free 
product of states, the normalization of square root type works in both cases since the 
number of summands in S(n) which give a non-zero contribution to the limits is in both 
cases of order n. 

Of course, the uniform boundedness assumption is satisfied if we take variables of 
type Xij{n) = Xij/^/n, where (Xij) is an infinite array of random variables whose 
distributions in the states (j)ij{n), respectively, are identical and do not depend on n. 
Although it is convenient to think of the Xjj(n) as if they were of this form, we want 
to study similar variables, whose variances are of order 1/n and stay the same within 
blocks whose sizes become infinite as n ^ go. 
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If the tuple {{ii,ji),---,{im,jm)) ^ {I x I)"^, where / is an index set, defines a 
partition tt = {tti, . . . , TTfe} of the set [m], i.e. {ip,jp) = {iq,jq) if and only if there exists 
r such that p,q e TTr, we will write 

Of course, if tt is a non-crossing pair partition, then each tt^ is a two-element set. We 
will also adopt the convention that if m is odd, then A/'C^ = and the summation 
over TT G A/'C^ gives zero. This allows us to state results for moments of the S{n) of all 
orders without distinguishing even and odd moments. 

Lemma 6.1. Let {An,4>n) be a sequence noncommutative probability spaces, each with 
an array {4>i,j{n)) defined as above and such that 

(1) for fixed n e N, the array (Xj .,(n)) is matricially free with respect to {(j)i^j{n)), 

(2) ^ij(n)(X,j(n)) = and (t)ij{n){Xlj{n)) = i\,.j{n) for any i,j and n e N, 

(3) the moments of the Xjj(n) are uniformly bounded with respect to (0jj(n)). 

Then 

^^J4S^{n))= J] v{7r,n) + 

where v{7r,n) = tr(y(7r, n)) is given by Definition 5.1 and corresponds to the variance 
matrix V{n) = {vij{nj). 

Proof. We have 

= 2 2 M^hJiin) ■ ■ ■ Xi^jjn)) 

-P(('l.Jl).---.(«iTi j'm ))=•"■ 

where Vm denotes the set of all partitions of [m]. Standard arguments allow us to 
conclude that for large n only non-crossing pair partitions give relevant contributions, 
namely 

(1) the moments of the Xij{n) satisfy the singleton condition with respect to the 
V'n.fc and thus with respect to V'n, namely 

^„(X,,,,,(n)...X,„,,„(n))=0 

if TT = P((^l, ji), . . . , {im,jm)) coutaius a singleton, 

(2) if TT has no singletons and contains blocks consisting of more than two elements, 
then 

{Xi,,j, (n) . . . Xi^j^ {n)) = 

(3) if m is even and tt is a crossing pair partition of [m], then the contribution from 
the corresponding mixed moments of the X^j is zero. 
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The above arguments are similar to those that would hold if random variables Xjj(n), 
where e I x I, were free. (1) follows from the definition of the matricially free 
product of states and the mean zero assumption. (2) holds since moments are uniformly 
bounded and the usual combinatorial argument (there are fewer than m/2 independent 
indices to sum over) can be used. In particular, this and (1) imply that if m is odd, 
then the contribution to the limit is of order 0{l/^/n). (3) is a consequence of the mean 
zero assumption and the definition of the matricial freeness (the mixed moments in our 
case vanish always when the corresponding moments vanish in the free case). 
Therefore, 

Now, suppose that m is even, tt g A/'C^ and the sequence of pairs ((ii,ii), . . . {im,jm)) 
is 'compatible with the matricial multiplication', i.e. such that if {ik,jk) and {ir,jr) 
label an inner-outer pair of blocks, then jk = v. If ttj = {r, r + 1} is a block which 
has no inner blocks, then we can 'pull out' the variance corresponding to that block, 
namely: 

i^n{Xi,,j,{n) . . . Xi^j^{n)) 

= Vi,,jMXi,j,{n) . . .^v_i,>_i(n)Xi^^2,>+2W • ■■Xi^jmin)) 
(mean zero assumption is used again). Continuing this procedure with other blocks 
which have no inner blocks, and summing over indices ii, ji, ■ ■ ■ ,im, jm, for which it 
holds that P{{ii,ji), . . . , {im,jm)) = we arrive at 

^ 2 (^) • • • ^k^Mm) ('^) + O 

where a{r) = if tt^ has no outer blocks (/cq labels the conditional block) and (7(r) = j 
if the nearest outer block of tt^ is labelled by j. Let us observe that we included in 
the above sum all possible labellings of the blocks of tt. This is done for convenience 
since it enables us to express the final result in terms of v(7r). More explicitly, we allow 
kQ,ki, . . . ,km to assume arbitrary values from the set [m] (in particular, they can all 
be equal), which produces certain terms which cannot be obtained from the summation 
over all ((ii, ji), . . . , (irmjm)) which define tt. For example, no nearest inner-outer pair 
of blocks can contribute vjjVjj, which appears in the above sum. However, all such 
terms are of order due to insufficient number of different summation indices (there 

are fewer than m/2 independent indices) and therefore they can be included in the sum 
without changing the asymptotics. Using Lemma 5.1, we obtain our assertion. ■ 

Lemma 6.2. Under the assumptions of Lemma 6.1 it holds that 

<j>n{S^{n))= 2 v^{'K,n) + 

where t'o(7r, n) = Tr(Vo(7r, n)) is given by Definition 5.2 and corresponds to the variance 
matrix V{n) = {vij{n)). 
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Proof. The proof is similar to that of Lemma 6.1. I 

In order to ensure existence of the limits of v{Tv,n) and vo(7r,n) as n ^ oo, we need 
to make some assumptions on the sequence of matrices {V{n))nm- For that purpose, 
introduce a matrix, called the dimension matrix 

D = dia.g{di,d2,. . . ,dr) e Dr{R) with Tr(D) = 1, 

where di,d2, . . . ,dr are positive real numbers. For given natural n, we associate with 
D the partition [n] = Ni u N2 u . . . u Nr, where 

Ni = [1, rii], N2 = [1 + ni,ni + 712], . . . , Nj. = [1 + rii + . . . + n^-i, n] 

and 

/ k \ /fc-l 
rik = E i^diTi] - E i^diU 

for each k e [r], with E{x) denoting the largest integer smaller or equal to x. Of course, 
the Nk are disjoint and naturally ordered intervals containing Uk natural numbers, 
respectively. Note that in the limit n ^ 00 we obtain rik/n d^ for each k. 
Assume now that each matrix V{n) has the block-form 

/ ^1,1 (n) ... Vi,rin) \ 
(6.4) V{n)= : ■■. : 

\ Vr,l{n) ... Vr,r{n) J 

where each block Vij{n) consists of the same number Uij/n. In other words, each V{n) 
is obtained from a real-valued matrix 



U = (uij) e Mr 

by repeating nj x Uj times each entry Uij at all entries of block Vij{n) and dividing it by 
n. Unless stated otherwise, we assume that Ujj > for any j and that Uij ^ for any 
i ^ j. The whole sequence {V{n))neN is built in such a way that proportions between 
sizes of these blocks are similar for all n and expressed in terms of the dimension matrix 
D (asymptotically, these proportions are given by the proportions between numbers dj). 

Assuming that the variance matrices V{n) are of the above block form, we can now 
state the standard and tracial central limit theorems, with limit distributions described 
in terms of traces of Section 5. 

Theorem 6.1. Under the assumptions of Lemma 6.1, if V{n) is of the block form (6.4) 
for each n eN, then 

(6.5) lim^45"(n))= V 6(7r), 

n—>oo '—^ 

for any m G N, where b{7r) = Tr(B(7r)D) and B(7r) is the diagonal matrix of Definition 
5.1 corresponding to tt and the matrix B = DU . 

Proof. Clearly, if m is odd, we get zeros on both sides of the above formula (we use 
our convention that in this case HC'^ = 0). The proof for m = 2k, where k e N, 
is based on Lemmas 5.1 and 6.1. If, in the combinatorial expression for v{7r,n), we 
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substitute for the matrix elements of V{n) the assumed block form, then, using the 
partition of the set of colors [n] = Ni^ . . .'u Nr, we can perform summations over 
the colorings which belong to each interval Nj separately. Thus, the contributions of 
various [n]-colorings of tt to the limit laws reduce to those corresponding to [r]-colorings 
and are described in terms of numbers Uj{7ii,f ), where i e [k], j e [r] and / g Fr{n) 
(the number j is the color of the conditional block). We have 

lim^;(7r,n) = lim ( Y! % V nf^x)Uj{'KiJ) . . .nf^k)Uj{i^kJ) 

\ je[r] /eF.(7r) 

je[r] /eF,.(7r) 

where tt -B(7r) is the matrix- valued function which corresponds to the matrix B = DU 
in accordance with Definition 5.1. In terms of matrix multiplication, the expression on 
the right hand side is obtained from that of Lemma 5.1 corresponding to matrix U 
by multiplying U from the left by the dimension matrix D and multiplying the whole 
product of matrices from the right by D. This proves our assertion. ■ 



Theorem 6.2. Under the assumptions of Lemma 6.2, if V{n) is of the block form (6.4) 
for each n eN, then 

(6.6) lim cPr^iS^'in)) = V bo{n), 

n^oo 

7reA/-C^ 

for any m e N, where n bQin) is the real-valued function of Definition 5.2 corre- 
sponding to the matrix B = DU. 

Proof. The proof is similar to that of Theorem 6.1 and is based on Lemmas 5.1 
and 6.2. The only difference is that we do not use conditional blocks to describe the 
colorings of all blocks of tt. Thus, the contributions of various [n]-colorings of n to the 
limit laws reduce to those corresponding to [r]-colorings and are described in terms of 
numbers u{7ri,f), where i e [k] and / g Fr{n). Namely, we have 

lim vo{7T, n) = lim ( ^ V nf,u{7ri, /) . . . nf^u{nk, f) ] 

n— >cx) n— >oo \ fl'^ I 

= Yi df,u{TTi,f)...df^u{7rkJ) = Tr(Bo(7r)) 

as n ^ CO, where tt Bo{7r) is the function defined by Definition 5.2, which proves our 
assertion. ■ 

7. Matricial semicircle distributions 

The results of Section 6 lead to combinatorial formulas for the asymptotic moments 
in the corresponding central limit theorems. In this Section we are going to express the 
limits in terms of their Cauchy transforms represented in the form of continued fractions. 
They play the role of the (standard and tracial) 'matricial semicircle distributions'. 
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For that purpose let us recall definitions of certain convolutions of distributions, or 
more generally, of probability measures. If is the reciprocal Cauchy transform of some 
probability measure fj, e Mr, then the K-transform of /i is given by K^{z) = z — Ff^{z). 
The boolean additive convolution fi u can be defined by the equation 

K^^^{z) = K^{z) + K^{z) 

where /U, z/ e Ai^ and z e C"*", respectively. In fact, this equation shows that the 
K-transform is the boolean analog of the logarithm of the Fourier transform [19]. 

We will also need another convolution, which reminds the monotone convolution [17], 
called the orthogonal additive convolution and defined by the equation 

K,U^) = K,{FAz)) 

where fj,,^ e A^r and z e C+. It was introduced in [12], where we showed that the 
above formula defines a unique probability measure on the real line. Moreover, if /i and 
v are compactly supported, both /i \- and /iwu are compactly supported. 

Using these convolutions, we will now define certain important continued fractions 
which converge uniformly on the compact subsets of C"*" to the K-transforms of some 
probability measures //jj e M.r- 

Lemma 7.1. For given B e Mr(R) with nonnegative entries, continued fractions of the 
form 



2jfe 



2-ip ' 



z — . . . 

where i,j e [r], converge uniformly on the compact subsets o/C"*" to the K-transforms 
of some iiij e Ai-R. with compact supports. 

Proof. Let us define a sequence of functions which approximate the Kij. Namely, set 
K^j'{z) = bij/z for any i,j e [r], which are the K-transforms of probability measures 
on R for any i,j (Bernoulli measures if bij > and if = 0). In order to use an 
inductive argument, let us establish the recurrence 

Kffiz) = 

for m ^ 1. If the are the K-transforms of some //p^ e M.r for any i and p, 

respectively, then the sums are the K-transforms of some 

and next, the K^ff' are the K-transforms of some 

where the the Bernoulli measures with K-transforms Kij(z) = bij/z, respec- 

tively. It is easy to see that all these measures are compactly supported. Moreover, the 
properties of the orthogonal additive convolution (Corollary 5.3 in [10]) say that the 
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moments of /Lt^"-' of orders ^ 2m agree with the corresponding moments of ^ for any 
n > m and any given Therefore, we have weak convergence 

w - lim n^^^ = fXij 

TO— >CX) '•' 

to some e M.^ for any These measures are also compactly supported since 
sup - ^ bi, j is finite. In turn, this implies that the corresponding Cauchy transforms (and 
thus K-transforms) converge uniformly to the Cauchy transform (K-transforms) of //.^j 
on compact subsets of C"*". This completes the proof. ■ 

We are ready to state a theorem, which can be viewed as the tracial central limit 
theorem for matricially free random variables. 

Theorem 7.1. Under the assumptions of Theorem 6.1, the 11)11- distributions of Sn con- 
verge weakly to the distribution given by the convex linear combination 

r 

where jij = /iij ty iJ,2,j ty . . . ty for each j = 1, . . . ,r and /lij is the distribution defined 
by Kij for any 

Proof. By Theorem 6.1, we have combinatorial formulas for the moments of the 
limit law in the tracial central limit theorem. The associated distribution extends to a 
unique compactly supported probability measure // on the real line since its moments 
are bounded by the moments of the Wigner semicircle distribution da with variance 
a = supi jbi j. Using the mulitplicative formula (5.4) for B{tt), we can formally write 
the Cauchy transform of jj, in the form 



k=0 



Tr {{z - K{z))-^D) 



which can be called the 'trace formula' for (7^, where 



00 



(7.1) K{z) = Y, 2 B{7r)z-"'^' 

k=l -neNCel^ 

is a diagonal-matrix-valued formal power series. Moreover, we will show below that 
each function Kj on its diagonal is, in fact, the K-transform of some e AIr. Then, 
the formal power series given by the trace formula is the Cauchy transform of as a 
convex linear combination of Cauchy transforms of probability measures. In fact, using 
the definition of i?(7r) and (5.3), we obtain the equation 

(7.2) K{z)=T{{z-K{z))-'B) 
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where K{z) = diag(i^i(2;), . . . , Kr{z)). By analogy with the scalar-valued case, we can 
find its solution in the form of a continued fraction. Namely, observe that each Kj{z) 
has the form of a formal Laurent series 

z z^ 

for some c_i,c_3, . . ., and therefore the above vector equation can be solved by succe- 
sive approximations. Namely, we set /Uj to be the (compactly supported) probability 
measure associated with the K-transform Kj{z) = ^^^^ j e [r], where 

the Kij are given by Lemma 7.1 for B = DU. These K-transforms solve (7.2). This, 
together with the trace formula for Gfj,, gives 

r 

i=i 

where G^.{z) = l/{z — Kj{z)) is the Cauchy transform of /Xj for j = 1, . . . ,r. That 
completes the proof. ■ 

Remark 7.1. The Cauchy transform of each /ij can be written as a continued fraction 
of the form 

G.A^) = 

^-Z^i 7 



^-Sfc 7 

which converges on the compact subsets of C"*". 

Next, we state a theorem, which plays the role of the standard central limit theorem 
for matricially free random variables. 

Theorem 7.2. Under the, assumptions of Theorem 6.2, the distributions of Sn con- 
verge weakly to the distribution 

1^0 = W 1^2,2 W . . . W llr,r 

where fijj is the distribution defined by Kjj for each j. 

Proof. By Theorem 6.2, we have combinatorial expressions for the limit moments 
M^. The proof is similar to that of Theorem 7.1 and is based on the trace formula for 
the K-transform of fj,o 

K,,{z)=Tr{{z-K{z))-'B,) 

derived from the definition of the function bo, which leads to the equation for the Cauchy 
transform 

^''^"^^ z-T.-K,,izy 

which completes the proof. ■ 
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Remark 7.2. The Cauchy transform can be written as a continued fraction of the 
form 

G,,{z) = — 

z - 2 



'•^ h- ■ 



z — . . . 

which gives a matricial extension of the continued fraction of the Wigner semicircle 
distribution. 

8. Decompositions in terms of subordinations 

In the one-dimensional case the limit distributions /iq and ^ are related by Proposition 
5.1. In particular, if each variance matrix V{n) has identical entries equal to one, both 
central limit theorems (standard and tracial) give the Wigner semicircle distribution 
with variance 1 (of course, the standard case also follows from free probability, whereas 
the tracial case is related to random matrices). 

In this Section we will analyze in more detail the limit distributions for the two- 
dimensional case, namely when each variance matrix V{n) consists of four blocks. They 
will be expressed in terms of two-dimensional arrays of distributions. Finding simple 
analytic formulas for the corresponding four-parameter Cauchy transforms and densities 
does not seem possible in the general case. However, we shall derive decomposition 
formulas for those measures in terms of s-free additive convolutions [12], which gives 
some insight into their structure (see also [18] for recent results on the multivariate case). 
The s-free additive convolution refers to the subordination property for free additive 
convolution, discovered by Voiculescu [24] and generalized by Biane [4]. As shown in 
[12] and [13], there is a notion of independence, called freeness with subordination, or 
simply s-freeness, associated with the s-free additive convolution and its multiplicative 
counterpart. 

Recall that the s-free additve convolution of /i, u e A4r is the unique probability 
measure □ e M.r defined by the subordination equation 

U S fjL = U V> {fj, [E I'), 

where o denotes the monotone additive convolution [17]. Equivalently, the above sub- 
ordination property can be written in terms of Cauchy transforms or their reciprocals. 

Using s-free additive convolutions and the boolean convolution, we obtain a decom- 
position of the free additive convolution of the form 

/ufflz/ = (/uffiz/) w (z/ffi/i), 

which allows us to interpret both s-free additive convolutions appearing here as (in 
general, non-symmetric) halves of fi^i/. We find it interesting that the limit distribu- 
tions in the two-dimensional case will turn out to be deformations of the free additive 
convolution of semicircle laws implemented by this decomposition. In other words, the 
subordination property and the associated convolutions give a natural framework for 
studying matricial generalizations of the semicircle law. 
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For simplicity, it will be convenient to use the indices-free notation for the two- 
dimensional matrix of K-transforms: 

a{z) b{z) \ ^ f K,,,{z) K,,2{z) 
c{z) d{z) J [k^Az) K,,,{z) 

and 



B 



where the square root is interpreted entry-wise. 

Moreover, we will distinguish two laws by special notations: we denote by (Tq the 
Wigner semicircle distribution with the Cauchy transform 

where the branch of \l z^ — 4q;^ is chosen so that \l z^ — ^o? > if z G R and z G (2q;, go), 
and by we denote the Bernoulli law with the Cauchy transform 

(^) = 

^ z — Y/z 

i.e. 0-,, = 1/2 + (5^). 

Finally, we will also use the boolean compressions of /i G A^r, where t ^ 0, defined by 
multiplying its K-transform by t, namely we define Tf/j, to be the (unique) probability 
measure on R, for which 

These transformations were introduced and studied in [7] and called 't-transformations' 
of /i. We allow t = 0, in which case T^ji = 6q. In particular, we shall use two-parameter 
boolean compressions of semicircle distributions, cr^,/? = TtCJa for t = {/3/aY, with 

being their Cauchy transforms, where the branch of the square root is the same as in 
the case of Ga^- 

Theorem 8.1. If a, l3,^f, 6 0, then the diagonal measures fijj defined by Kj j, where 
1 ^ j ^ 2, have the form 

with the non-diagonal measures given by fii^2 = ^tA*i,i o-nd //2,i = 7"sA*2,2; where t = 
{p/af and s = {^/5f, 

Proof. It is easy to see that the following algebraic relations hold: 

a{z) = — and d{z) 



z — a{z) — c{z) z — b{z) — d{z) ' 

P^ 7^ 
b{z) = -— and c(z) 



z — a{z) — c{z) z — b{z) — d{z) 
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Thus, b{z) = ta{z) and c{z) = sd{z), which gives fj,i^2 = TffXi^i and //2,i = 2^sA*2,2- In 
turn, from the equation for a{z), we get 

and thus /xi^i = cTq, |— ^2,1- In a similar manner we obtain 112,2 = crs \- /^i,2- Therefore, 
we arrive at the equations 

A*2,2 = (TSV- {Ttaa h Tsll2,2) 

since Tt(/x h i^) = (TtA*) l~ ^- In order to express the fXjj in terms of s-free additive 
convolutions, we need to use the properties of the orthogonal convolution. We have 
shown in [12] that the moment of order k oi fi \- depends on the moments of orders 
^ A; of /i and the moments of orders ^ k — 2 oi u. This leads to the conclusion that for 
any compactly supported 11, u e M.^ and the associated sequence of measures {(i hm i^), 
defined recursively by 

A* l-TO = I- ('^ l-m-i y") with hi = /U h i^, 

we have weak convergence w — lim^^oo n hm = □ z^. If we take ji = T^Ua and 
u = TgOJs (these measures are compactly supported), we get the desired formulas. ■ 



Corollary 8.1. The measures /j^o, fii, IJ-2 can be decomposed as 

Hi = Ti/t{aa,(s ffl o-s,^) w {'Js,^ ffl o-a,/?) 

where the assumptions and notations are the same as in Theorem 8.1. 

Proof. These decompositions follow immediately from Theorems 7.1, 7.2 and 8.1. ■ 

Remark 8.1. The formulas for the diagonal measures fij^j in the proof of Theorem 8.1 
remind those for 2-periodic continued fractions if we take t = s = 1. The latter are of 
the same form, except that the semicircle distributions are replaced by much simpler 
Bernoulli laws. Nevertheless, \ft = s = 1, the formulas for the jij take a simple form 

where j = 0,1,2. By Theorem 7.1 and Corollary 8.1, the same formula holds for yu. 
Therefore, all measures /io, /ii, /U2, fi can be viewed as deformations of the free additive 
convolution of two semicircle distributions, implemented by means of boolean compres- 
sions. 

Let us consider now the situation in which some of the numbers a, /9, 7, 5 vanish. 
Suitable formulas can be derived algebraically, as we did in the proof of Theorem 8.1. 
However, one can also obtain the same results by taking weak limits in the formulas 
for the measures yUjj, using the fact that all measures involved have compact supports. 
For that purpose, let us state a few useful facts about weak limits which will be of 
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interest to us. Then, we consider eight cases, to which the remaining cases are similar 
(for instance, a = /3 = is similar to 7 = (5 = 0). 

Proposition 8.1. Let t = P'^/a'^, where a, /3 > 0, and let /i e Aim be compactly 
supported. 

(1) Ifa^ 0+, then 

(a) w — lim cTq = Sq, 

(b) w - lim (aa,^) = Kg, 

(c) w - lim {Ti/t{aa,(3 ffl A*)) = <^o- 

(2) If 13^0+ , then 

(a) w - lim cr«,/3 = So, 

(b) w - lim (Ttfi) = So, 

(c) w - lim (Ti/t(c7«,/3 ffl /u)) = (Tc, h /u. 

Proof, li a ^ 0"*", then K^^(z) 0, which proves (la). Here, as well as in the 
remaining cases, convergence is uniform on compact subsets of C"*". Moreover, K^^ ^ = 
(3^/{z - a^K„^) = K^^, which gives (lb). In turn, if /9 ^ 0+, then K^^ ^{z) = 

P'^/{z — Ka^(z)) 0, thus also w — lim aa.a = So, which proves (2a). If, in addition 
t ^ 0"*", then Tffx So for any fj, e A4k, which proves (2b). Finally, 

ri/t((7a,/3ffl/u) = Ti/t{Ttaa h {lJ,[Ecra,f})) = CTah (// ffl (Ta,;?) 

and thus the right hand side tends weakly to So as a ^ 0^, which gives (Ic), and tends 
weakly to cr^ h (^u □ So) = aa h ^i, which gives (2c). This holds for any fj, e Adm, and 
we also use the right unit property of S with respect to the s-free additive convolution, 
namely /i ffi = A*- ' 

Corollary 8.2. If some of the entries of the matrix A vanish, we can distinguish eight 
different cases, for which the distributions fiij are given by Table 1. 

Proof. If ajj = 0, then yUjj is the Dirac delta, which easily follows from the algebraic 
equations for the corresponding K-transforms. An alternative proof can be given by 
taking weak limits of the formulas of Theorem 8.1, as we proceed with the remaining 
measures. Thus, if 5 ^ 0+, then 

Atl,l = W - lim Tyt{(Ta,/3 ffl (75,7) = Tiit{(Ja,p ffl K^) 

Ati,2 = w - lim((7«,;3 □ as,^) = aa,p ffl k^, 
1^2,1 = w - lim ((75,^ ffi (7«,/3) = n^m cra,p: 

by (lb) of Proposition 8.1, which proves the first case in Table 1. The remaining cases 
are proved in a similar manner. ■ 

In the case of arbitrary matrix A, finding the four-parameter densities of /io and 
is unwieldy. Below we shall just consider two special cases, in which we can find nice 
formulas for these measures for matrices A of arbitrary dimension. These two cases are 
of special interest since they are associated with (aymptotic) freeness and (asymptotic) 
monotone independence. 
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Table 1. Distributions /Ujj in the case A has zero entries 



ai,i 


ai,2 


^2,1 


^2,2 


1^1,1 




1^2,1 


fJ'2,2 


a 


P 


7 







(^a,f3 ffl /«7 




So 







7 





So 


Hp ffl S 




So 


a 





7 


6 













(3 





5 


So 




So 


as \- Kfs 








7 


5 


So 


5o 






a 








5 


aa 


^0 


So 




a 















So 


So 





/3 








So 




So 


So 



Proposition 8.2. If A is a square r -dimensional matrix with identical positive entries 
aj in the j-th row, then 

l^j = Cr^i S 0"a2 S • • • ffl "^ar 

for each j e [r], and the jij coincide with fi and /lo- 

Proof. Since the columns of A are identical, the functions Kij are the same for all 
j's. Denote them Li = Kij, where i,j e [r]. Moreover, 

for any i e [r]. Therefore, 21=1 ^ii^) ^he K-transform of the measure 

(^ai ffl (^aa ffl ■ ■ ■ ffl CTar 

and since Ki^.{z) = XiLi ^hji'^) ^ SLi -^«(-^)' ^^e proof for /ij is completed. It is then 
easy to see that we get the same result for /i and /iq. ■ 

Proposition 8.3. If A is a lower-triangular r -dimensional matrix with identical positive 
entries aj in the j-th row below and on the main diagonal, then 

for each j e [r] . Moreover, /iq = /ii and fj, is the convex linear combination of the /ij as 
in Theorem 7.1. 

Proof. As in the proof of Proposition 8.2, note that the /Ujj do not depend on j and 
thus we can set Li = Kij for any i ^ j. If i = r, we have 

^'■^^^ = z-i{zy 

using the continued fraction for the Kj.j of Lemma 7.1. Therefore, firj = <^ar for 
j ^ r. Next, we have 

which leads to 

r 

L^{z) = K,^^{z - U{z)) 
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which gives the orthogonal decomposition of fj,k,j, 
for any j ^ k < r. Now, we claim that 

for any 1 ^ j ^ i ^ r. Clearly, it holds for z = r and any j ^ r since we have already 
shown that = (Ta^ for any j ^ r. Suppose now that this formula holds for z > and 
any j ^ i. We will show that it holds ioi i = k and any j ^ k. Using the orthogonal 
decomposition of //^j given above and the inductive assumption, we obtain 

= (Ta^ h {cTa^^, t> (. . . t> a« J . . .) ■ 

However, for any /U, e M.r, we have a simple relation 

(/U|-i^)tyi^ = /uc>z^ 

which gives 

and the desired expression for fXj. In a similar manner we obtain /xq and fi. ■ 

Example 8.1. If a = S ^ and /? = 7 = 0, then we can use Table 1 to obtain 
jjio = Oa ^ (Ta, which is the arcsine law with Cauchy transform ^^^0(2;) = l/\fz^^^^o? 
whereas jJL is the Wigner semicircle distribution ctq. In turn, if a = 5 = and /3 = 7 7^ 0, 
then jiQ = whereas /j, = crp. 

9. Weighted binary trees and Catalan paths 

In thie Section we show how to express the limit distributions in terms of walks on 
weighted binary trees, or equivalently, in terms of weighted Catalan paths. The binary 
tree serves here as an example of the strongly matricially free Fock space. 

The usual framework which gives a description of distributions in terms of walks on 
graphs is the following. Let W{n) denote the set of root-to-root walks of lenght n on a 
rooted graph {Q,e) and let fi be the spectral distribution of {G,e), i.e. the distribution 
given by the moments of the adjacency matrix A{Q) in the state (p associated with the 
vector 6e on the space of square intergrable functions on the set V{Q) of the vertices 
of Then the n-th moment of A{Q) in the state (p is equal to the cardinality of the 
set W{n). In particular, it is well known that the moments of the Wigner semicircle 
distribution of variance 1 can be expressed in terms of walks on the half-line (Ti,e) 
with the first vertex denoted by e and chosen as the root. 

For many distributions we have to use a more general framework, in which the mo- 
ments of these distributions are expressed in terms of root-to-root (random or, more 
generally, weighted) walks on some rooted graph, except that to each walk w on this 
graph we have to assign a real- valued weight ^{w). Then we can write 

'weW(n) 

for any n ^ 1, where fj, is the considered distribution. In particular, we obtain the 
moments of Ua for any a > by putting ^{w) = a^, where n = \w\ is the lenght of w. 
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Figure 3. Binary tree with a matricial weight function 

In the cases which are of interest to us, the weight function is first defined on the 
set of edges E{Q) of Q and then is extended to W = Un^i^l*^) multiplicativity. 
Namely, if we are given a mapping ^ : E[Q) R, we set 

eH = e(i?iX(i?2)...e(K) 

where w = [Ei, E2, . . ■ , E^) and Ei, E2, ■ ■ ■ , En are the edges of w (we choose to describe 
walks on graphs as sequences of edges). Such extension, by abuse of notation denoted 
also by ^, will be called multiplicative. 

For instance, it is easy to see that the moments oi /i = a^So's can be expressed in 
this form. It is enough to take the free product of two half-lines, which is the binary 
tree (T2,e) with root e. Let us color this graph in the natural way, namely each edge 
which belongs to a copy of the first half-line is colored by 1 and each edge which belongs 
to a copy of the second half-line is colored by 2. Then the above formula holds for the 
moments of fj, if we take ^ (E) = a whenever E is colored by 1 and ^ (E) = 5 whenever 
E is colored by 2. 

We will demonstrate below that we can express our distributions Ho, Hi, H2 in a similar 
form (in particular, we can use the binary tree), except that the weight function ^ will 
depend on all four parameters which appear in the matrix A. Before formulating the 
theorem, let us introduce special weight functions related to matricial freeness. 

Definition 9.1. Let be a r-ary rooted tree with root e and let A e M,.(M). The 
weight function ^ : E{Tp) R which assigns the entries of A to the edges of is 
called matricial if, for any pair of edges Ei,E2 e E{Tr), incident on the same vertex 
and such that Ei is the 'father' of E2, the following implication holds: 

^(£■1) = ttij for some i,j =^ ^(-^2) = 0'k,i for some k. 

The unique multiplicative extension of this weight function to the set of all walks on 
Tr will also be called matricial. 

We specialize to p = 2 and the binary tree. Note that any matricial weight function ^ 
on the binary tree is uniquely determined (up to equivalence) by the set of those entries 
of the matrix A which are assigned to the set {Ei,E2} of two edges incident on the 
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root of tree, called the initial weights. In order to establish a connection with our limit 
distributions, we will assume, as in Section 8, that A is the 'square root' of B = DU 
of Section 6. Then, in particular, the binary tree with the matricial weight function 
associated with A and the initial weights {a, 5}, shown in Figure 3, describes jiQ as we 
show below. 

Finally, recall after [1,12] that if (^i,ei) and (^2,62) are two locally finite simple 
graphs and /xi and 112 are the asscociated spectral distributions, then the s-free product 
of Qi and Q2 (in that order) can be interpreted as this half of the free product Qi * Q2 
which 'begins' (starting from the root) with a copy Qi. Moreover, the associated spectral 
distribution is given by /xiffi/X2 Let us add that a similar result holds in the multiplicative 
case: both the s-free multiplicative convolution and the associated s-free loop product 
of graphs were introduced in [13]. 

Below we shall use the s-free product of half-lines, which are (left and right) halves 
of the binary tree. 

Theorem 9.1. Let ^q, ^1,^2 be the multiplicative matricial weight functions on the set 
of walks on T2 associated with matrix A and the initial weights {a, S}, {(3, S} and {a, 7}, 
respectively. Then 

(9.1) M,^{n)= J] ^j{w) 

'weW(n) 

for j = 0, 1, 2 and any n e N, where W{n) denotes the set of root-to-root walks on T2 
of lenght n. 

Proof. We need to translate the result of Corollary 8.1 to the language of graphs. 
It is well-known that the moments of (Tq can be interpreted in terms of walks on the 
half-line with the weight a assigned to each edge. The boolean compression Tt of cTq 
changes only the weight assigned to the edge incident on the root, namely it multiplies 
it by = P/ce, which can be illustrated as 

a a a a a p a a a a 
m • • • • ~^ • • • • • 

Now, the s-free additive convolutions of compressed semicircle distributions which ap- 
pear in the decompositions of Corollary 8.1, namely 

(7a,p\Bcr 5,-y and (Ts,'^\E(Ta,p 

are the spectral distributions of the s-free products of these half-lines (taken in two 
different orders) . These turn out to be the spectral distributions of the two halves of 
the binary tree T2 with the weight function defined by the initial set 7}. In order to 
obtain /io, we still need to apply Tyt and Ti/g, respectively, to the left and right halves 
of the tree, which amounts to changing the initial weights from P and 7, respectively, 
to a and 5. As a result, we obtain the weight function associated with A and the initial 
set {a, 5}. This proves the statement concerning //q (the corresponding weight function 
on the binary tree is shown in Figure 3). The cases of /Ui and ^2 are very similar (at 
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Figure 4. A weighted Catalan path associated with A 

the end, a boolean compression is applied to only one half of the binary tree). ■ 

Another geometric interpretation of Corollary 8.1 can be given in terms of Catalan 
paths. In order to define a Catalan path, we begin with a function / : [2n] [n], such 
that /(O) = /(2n) and — /(i — 1)| = 1 for any 1 ^ z ^ 2n, and then we define a 
Catalan path as its unique extension / : [0, 2n] [0, n] (by abuse of notation, denoted 
by the same symbol) obtained by connecting each (i — l,f{i — 1)) with {i,f{i)) with 
a segment, where 1 ^ i ^ 2n. Clearly, each Catalan path consists of segments of two 
types: 'rises' Ri, R2, ■ ■ ■ , Rn, and 'falls' Fi, F2, . . . , F^. Moreover, to each 'rise' Rj there 
corresponds the closest 'fall' F^(j) lying to the right of Rj and on the same (vertical) 
level. 

There is a natural mapping from the set W{2n) of walks of lenght 2n on the binary 
tree and the set C(n) of Catalan paths of lenght 2n. In order to rephrase Theorem 8.1 
in terms of Catalan paths, we need to take the sets of weighted Catalan paths, by which 
we understand pairs (f,^), where / is a Catalan path and ^ is a real- valued weight 
function defined on the set of segments of /. The multiplicative formula 

^if)=C{Rl)...mn)^{Fl)...^{Fn) 

assigns the corresponding weight to /. Weighted Catalan paths of special type defined 
below allow us to rephrase Theorem 8.1. 

Definition 9.2. A weighted Catalan path {f,C) of lenght 2n is called matricial if ^ 
assigns entries of A G Mp(M) to the segments of / in such a way that the following 
implications holds: 

C{Ri) = ttij for some i,j =^ ^(-^2) = o,k,i for some k, 

CiFi) = ttij for some i,j =^ ^(-^2) = for some k, 

for any two consecutive 'rises' Ri- R-i and two consecutive 'falls' ^1,-^2, and the same 
weights are assigned to Ri and for each i e [n\. 

We will consider below the set of matricially weighted Catalan paths of lenght 2n 
associated with the matrix A. In order to rephrase Theorem 9.1, using weighted Catalan 
paths, we need to restrict the set of weights which can be assigned to the first segment 
of each path (by analogy with trees, we call them initial weights). An example of a 
weighted Catalan path contributing to /Uq is given in Figure 4. 
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Corollary 9.1. Let Co{n), Ci{n) and 6*2(71) be the sets of matricially weighted Catalan 
paths associated with A, with the initial weights {ck,S}, {P,5} and {7,0;}, respectively. 
Then 

M,^(2n)= J] ^(/) 

(/,OeC,(n) 

for j = 0, 1, 2 and any n eN. 

Proof. This is an easy consequence of Theorem 9T since there is a natural bijecton 
between each Cj{n) and the pair {W{2n),^j) for any n. ■ 
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